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We give a review of recent work aimed at understanding the dynamics of gravitational collapse in 
quantum gravity. Its goal is to provide a non-perturbative computational framework for understand- 
ing the emergence of the semi-classical approximation and Hawking radiation. The model studied 
is the gravity-scalar field theory in spherical symmetry. A quantization of this theory is given in 
which operators corresponding to null expansions and curvature are well defined. Together with the 
Hamiltonian, these operators allow one to follow the evolution of an initial matter-geometry state 
to a trapped configuration and beyond, in a singularity free and unitary setting. 
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CN \ I. INTRODUCTION 

c. 

One of the outstanding problems in theoretical physics is the incomplete understanding at the quantum level of the 
' ' formation, and subsequent evolution of black holes in a quantum theory of gravity. Although a subject of study for 
00 over three decades, it is fair to say that, in spite of partial results in string theory and loop quantum gravity, there 
is no widely accepted answer to many of the puzzles of black hole physics. This is largely because there has been no 
O , study of quantum dynamical collapse in these approaches. Rather, progress has focused mainly on explanations of the 
O^' microscopic origin of the entropy of static black holes from state counting. A four-dimensional spacetime picture of 
^ ' black hole formation from matter collapse, and its subsequent evolution is not available in any approach to quantum 
gravity at the present time. 

This paper summarizes an attempt to address this problem in the context of Hawking's original derivation of black 
hole radiation: spherically symmetric gravity minimally coupled to a massless scalar field. This is a non-linear 2d 
J> ■ field theory describing the coupled system of the metric and scalar field degrees of freedom. Gravitational collapse 
\ in the classical theory in this model has been carefully studied numerically [ll, [2] , but its full quantization has never 
been addressed. 

' Hawking's semi-classical calculation |3| uses the eikonal approximation for the wave equation in a mildly dynamical 
" [ . background, where the dynamics centers on the surface of a star undergoing collapse. The essential content of it is 
the extraction of the phase of the ingoing mode from an outgoing solution of the scalar wave equation as a classically 
\ collapsing star crosses its Schwarzschild radius. According to this calculation, emitted particles appear to originate 
near the event horizon. This means that an emitted particle observed by a geodesic observer is transplankian at 
creation origin due to the gravitational redshift (which is infinite at the horizon). Its back reaction is therefore not 
negligible, bringing into question the entire approximation. 

It is likely that a complete understanding of quantum evolution in this system will resolve all the outstanding 
^ problems of black hole physics in the setting in which they originally arose. The following sections contain a summary 
' of the work described in refs. [1, H, H, Bi- 
ll. CLASSICAL THEORY 
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The phase space of the model is defined by prescribing a form of the gravitational phase space variables qab and 
tt"'', together with falloff conditions in r for these variables, and for the lapse and shift functions N and TV", such that 
the ADM 3-1-1 action for general relativity minimally coupled to a massless scalar field 



^"'qab + -NH- NXa (1) 
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is well defined. The constraints arising from varying the lapse and shift are 

+47rG (^-^P| + ^q'^'da^dbc^^ = 0, (2) 
Ca = D.irl - P^da^ ^ 0, (3) 

where and 7?. is the Ricci scalar of Qab- The falloff conditions imposed on the phase space variables are 

motivated by the Schwarzschild solution in Painleve-Gullstand (PG) coordinates, which itself is to be a solution in the 
prescribed class of spacetimes. These conditions give the following falloff for the gravitational phase space variables 
(for e > 0) 

Qab^eab+ ^3/,^^ +0[r ), TT = ^3/2 + +^(^ )' (4) 

where f°'^,g°'^, h'^'' are symmetric tensors, 7r°^ = n^-'^/y/q, and g = detga?,. 
In this general setting we use the parametrization 

Qab = A{r,t)'^ UaUb^ ^ {Cab - naUb) (5) 

for the 3-metric and conjugate momentum for a reduction to spherical symmetry, where Cab is the flat 3-metric and 
n° — x°'/r. Substituting these into the 3+1 ADM action for general relativity shows that the pairs {R,Pb) and 
(A, Pa) are canonically conjugate variables. We note for example the Poisson bracket 



(7) 



which is the bracket represented in the quantum theory (described below). 

The falloff conditions induced on these variables from ([U, together with those on the lapse and shift functions, 
ensure that the reduced action 

5i? ^ y" dtdr {PrR + PaA + P^cf) - constraint terms) + j dt{N'' KPk) (8) 

is well defined. This completes the definition of the classical theory. 

At this stage we perform a time gauge fixing using the condition A = 1 motivated by PG coordinates. This is 
second class with the Hamiltonian constraint, which therefore must be imposed strongly and solved for the conjugate 
momentum Pa. This gauge fixing eliminates the dynamical pair (A, Pa), fixes the lapse as a function of the shift, and 
leads to a system describing the dynamics of the variables {R, Pr) and {4>, P^) (B]. The reduced radial diffeomorphism 
generator 

Cred = P'a{R. 4>, Pr, n) + P^^' + PrR' = (9) 

remains as the only first class constraint. It also gives the gauge fixed Einstein evolution equations via Poisson 
brackets, for example 4> — {'/'j/ N^Cred}- 



III. QUANTUM GRAVITY 



The quantization route we follow is unconventional in that field momenta are not represented as self-adjoint op- 
erators; rather only exponentials of momenta are realized on the Hilbert space. This is similar to what happens 
in a lattice quantization, except that, as we see below, every quantum state represents a lattice sampling of field 
excitations, with all lattices allowed. This quantization allows definitions of bounded inverse configurations operators 
such as l/x, which for quantum gravity leads to the mechanism for curvature singularity resolution described below. 
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A quantum field is characterized by its excitations at a given set of points in space. The important difference from 
standard quantum field theory is that in the representation we use, such states are normalizable. A basis state is 

\^^T.k'^kPR(xk) ^^^L^T.,biP^(vi)) ^ \ai...aN,-M...bN2), (10) 

where the factors of L in the exponents reflect the length dimensions of the respective field variables, and 0^,6; are 
real numbers which represent the excitations of the scalar quantum fields R and (p at the radial locations {xk\ and 
{yi\. The inner product on this basis is 

(ai . . . aAT^; 5i, . . . hN2\a'i ■ ■ ■ a'j^^;b[ . . . b'^J = Sa^y^ . . . ^jva,!)^^ ' 

if the states contain the same number of sampled points, and is zero otherwise. 
The action of the basic operators are given by 

Rf |ai . . . aNi;bi . ..bj^^) = L'^'^akf{xk)\ai . . . aN^;bi . . . bN^)^ (H) 

k 

e^>^jPRi^j)\ai . . .aNi^bi . . .bN^) = |ai...,aj - Aj, . . . oatj; 61 . . . ^at^), (12) 

where aj is if the point Xj is not part of the original basis state. In this case the action creates a new excitation at 
the point Xj with value — Aj. These definitions give the commutator 

= -\f{x)L^e^>^^). (13) 

Comparing this with ([7]), and using the Poisson bracket commutator correspondence {,}<-> ih[ , ] gives 
where Ip is the Planck length. There are similar operator definitions for the canonical pair {(f), P^). 



Rj^etXPa(x) 



A. Singularity resolution 

To address the singularity avoidance issue, we first extend the manifold on which the fields R etc. live to include 
the point r = , which in the gauge fixed theory is the classical singularity. We then ask what classical phase space 
observables capture curvature information. For homogeneous cosmological models, a natural choice is the inverse scale 
factor a{t). For the present case, a guide is provided by the gauge fixed theory without matter where it is evident that 
it is the extrinsic curvature that diverges at r = 0, which is the Schwarzschild singularity. This suggests, in analogy 
with the inverse scale factor, that we consider the field variable 1/i? as a measure of curvature. A more natural choice 
would be a scalar constructed from the phase space variables by contraction of tensors. A simple possibility is 

The small r behaviour of the phase space variables ensures that any divergence in tt is due to the 1/i? factor. We 
therefore focus on this. A first observation is that the configuration variables R{r, t) and (j){r, t) defined at a single point 
do not have well defined operator realizations. Therefore we are forced to consider phase space functions integrated 
over (at least a part of) space. A functional such as 



Rf = / drfR (15) 
"'0 

for a test function / provides a measure of sphere size in our parametrization of the metric. We are interested in 
the reciprocal of this for a measure of curvature. Since i? ~ r asymptotically, the functions / must have the falloff 
/(r) ~ -r-2-e foj- Jij %Q be -^gll defined. Using this, it is straightforward to see that 1/i?/ diverges classically for small 
spheres: we can choose / > of the form / ^ 1 for r << 1, which for large r falls asymptotically to zero. Then 
Rf ^ and 1/i?/ diverges classically for small spheres. 

A question for the quantum theory is whether 1/i?/ can be represented densely on a Hilbert space as 1/i?/. This is 
possible only if the chosen representation is such that i?/ does not have a zero eigenvalue. If it does, we must represent 
1/i?/ as an operator more indirectly, using another classically equivalent function. Examples of such functions are 
provided by Poisson bracket identities such as 
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|i?/| V*G/(r 



g-^PR{r) 



y^,e^^-M|) , (16) 
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where the functions / do not have zeroes. The representation for the quantum theory described above is such that 
the operator corresponding to Rf has a zero eigenvalue. Therefore we represent 1/Rf using the r.h.s. of pB]) . The 
central question for singularity resolution is whether the corresponding operator is densely defined and bounded. This 
turns out to be the case. 

Using the expressions for the basic field operators, we can construct an operator corresponding to a classical 
singularity indicator: 



1 



\Rf\ Vlfi^i) 



-iPR{Xj) 



(17) 



The result is that basis states are eigenvectors of this operator, and all eigenvalues are bounded. This is illustrated 
with the state 



(18) 



which represents an excitation oq of the quantum field Rf at the point of the classical singularity: 



Rf\Sa„) = (2/|,)/(0)ao |S, 



(19) 



|Sao) 



(|aor/2-|ao-ir/2)'|S.„) 



which is clearly bounded. This shows that the singularity is resolved at the quantum level. In particular if there is 
no excitation of Rf at the classical singularity, ie. ao = 0, the upper bound on the eigenvalue of the inverse operator 
II. 



is 2/lj 



B. Quantum black holes 



The event horizon of a static or stationary black hole is a global spacetime concept. It does not provide a useful 
local determination of whether one is inside a black hole. The fundamental idea for defining a black hole locally is 
that of a trapped surface, first introduced by Penrose. One considers a closed spacelike 2-surface in a spacetime, 
and computes the expansions 9+ and 0_ of outgoing and ingoing null geodesies emanating orthogonally from the 
surface. If 0+ > and 0_ < 0, the surface is considered normal. On the other hand if 6*+ < and 0_ < 0, the 
surface is called trapped. This provides a criterion for subdividing a spacetime into trapped and normal regions. The 
outer boundary of a trapped region may be considered as the (dynamical) boundary of black hole, also known as 
the "apparent horizon" in numerical relativity. It is a function computed in classical numerical evolutions to test for 
black hole formation. Similarly, a setting for studying quantum collapse requires an operator realisation of the null 
expansion "observable," and a criterion to see if a given quantum state describes a "quantum black hole." 

The classical expansions in spherical symmetry are the phase space functions 

e'i = - ^ (2i?2 AA' ± Pa + iA^RR') . (20) 

Given phase space functions on a spatial hypersurface S, the marginal trapping horizon(s) are located by finding the 
solution coordinates r = (i = 1 • • • n) of the conditions ^-|- = and 0_ < 0, (since in general there may be more 
than one solution). The corresponding radii Ri = R{ri) are then computed. The size of the horizon on the slice E is 
the largest value in the set {Ri}. 

Since only translation operators are available in our quantization, we define Pa indirectly by 

where < A ^ 1 is an arbitrary but fixed parameter, and U\ denotes exp{iXPA/L). This is motivated by the 
corresponding classical expression, where the limit A ^ exists, and gives the classical function Pa. A is perhaps best 
understood as a ratio of two scales, A = Ip/loj where Iq is a system size. As for a lattice quantisation, it is evident 
that momentum in this quantisation can be given approximate meaning only for A <C 1. A is also the minimum value 
by which an excitation can be changed. 
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Definitions for the operators corresponding to R' and A' are obtained by implementing the idea of finite differencing. 
We use narrow Gaussian smearing functions with variance proportional to the Planck scale, peaked at coordinate points 
fk + dp, where < e <C 1 is a parameter designed to sample neighbouring points: 



fe{r,rk) = exp 



{r - Tk - dpY 
21% 



Denoting Rf^ by i?e for this class of test functions we define 



(22) 



R'{rk) := 7^ iRe-Ro). (23) 



Ipe 

Putting all these pieces together, we can construct the desired operators 

hirk) = -^^oAo (a. - Ao) T ^ {Ux - ul) Ao - -^^Ro (Rc - i?o) , (24) 

which have a well defined action on the basis states. 

In analogy with the classical case, we propose that a state \^!) represents a quantum black hole if 

(*|^+(rfe)|*) =0, and {^\Lirk)\^) < 0. (25) 

for some r^. The corresponding horizon size is given by Rh — {'^\R{rk)\'^). 

This definition is utilised as follows: Given a state with field excitations at a set of coordinate points {r,;}, one 
would plot the expectation values in Eqn. pS)) as functions of (^|i?(rfc)|\l/), and locate the zeroes, if any, of the 
resulting graph. The resulting "quantum horizon" location is invariant under radial diffeomorphisms because these 
act on states to shift the coordinate locations of field excitations, but leave the expectation values unchanged - the 
graph is a physical observable. 

It is straightforward to construct explicit examples of states satisfying these quantum trapping conditions. Some 
examples appear in [6]. The quantum horizons so determined are not sharp since {0'^) ^ 0. 



IV. SUMMARY AND OUTLOOK 



The results so far from this approach to understanding black hole formation in quantum gravity are threefold: (i) 
A quantization procedure which allows explicit calculations to be done, (ii) a test for black holes in a full quantum 
gravity setting which makes no use classical boundary conditions at event horizons, and (iii) singularity free and 
unitary evolution equations using the Hamiltonian defined in 

The main computational challenge is to use the formalism to explicitly compute the evolution of a given matter- 
geometry state until it satisfies the quantum black hole criteria, and then to continue to the evolution to see if and 
how Hawking radiation might arise. This work is in progress. 
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